Dissociation of composite vortex lattice in two-band superconductors in the flux-flow state 
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In multiband superconductors, each superconducting condensate supports vortices with fractional quantum 
flux. In the ground state, vortex in different bands is spatially bounded together to form a composite vortex, 
carrying one quantum flux O - Here we predict dissociation of the composite vortices lattice in the flux flow 
state due to the disparity of vortex viscosity and flux of the vortex in different bands. For a small driving current, 
the composite vortices starts to deform, but the constituting vortices in different bands move with the same 
velocity. For a large current, composite vortices dissociate and vortices in different bands move with different 
velocities. The dissociation transition shows up as increase of flux flow resistivity. In the dissociated phase, 
Shapiro steps are developped when an ac current is superimposed with a dc current. 
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Multiband superconductivity is realized in many supercon- 
ductors, such as the well studied V3Si, Nt^Se, the recently 
discovered MgB2 |1| and iron-based superconductors J3J- In 
these superconductors, electrons in different bands are cooled 
into distinct superconducting condensates, which interact with 
each other through interband tunneling. Multiband supercon- 
ductivity may also exist in the proposed liquid hydrogen under 
pressure where both the proton and electron bands contribute 
to superconductivity|3]|4]. In this case the interband tunneling 
is absent. The multiband superconductors attract considerable 
attention recently and phenomena unique to multiband su- 
perconductors, which do not have counterpart in single-band 
superconductors, have been predicted [ 5—1-4- 1 and observed 
experimentallv lflSHTTII . 

In multiband superconductors, there are several supercon- 
ducting gaps \Af,\ expO'6^), characterizing the quasiparticle ex- 
citation for superconducting condensate in each band respec- 
tively. Because the gap function is a complex function, each 
condensate thus supports vortex excitation with fractional 
quantized flux. [7 1 The energy of a single fractional quantized 
vortex diverges logarithmically due to the counter flows of dif- 
ferent condensates that have no charge transfer and are not 
coupled with magnetic fields. Thus it is thermodynamically 
unstable. Stable fraction quantized vortex is predicted to exist 
in mesoscopic superconductors! 18, 19 1 or near the surface of 
multiband superconductors . Il20l 

Vortices in different condensates appear simultaneously 
when an external field is applied to multiband superconduc- 
tors. There are both interband and intraband vortex interac- 
tion. Vortices with the same polarization in the same conden- 
sate repel each other through the exchange of massive pho- 
ton. The vortices in different condensates interact repulsively 
due to the magnetic interaction. Meanwhile they attract each 
other due to the coupling to the same gauge field. The latter 
is more important and the net interaction of vortices in dif- 
ferent condensates is attractive. They also attract with each 
other due to the interband tunneling in superconducting chan- 
nel. Therefore vortex in different condensates in the ground 
state is bounded and their normal cores are locked together 
to form a composite vortex with the standard integer quan- 



tum flux <t>o = hc/(2e). It is an interesting question whether 
the composite vortex can dissociate. The melting of com- 
posite vortex lattice due to thermal fluctuations was studied, 
and it was found that the vortex lattice in the superconducting 
condensate with weaker phase rigidity melts first as tempera- 
ture increases, while the lattice ordering in other condensates 
remains. IfTTII 

Here we consider the dissociation of composite vortex lat- 
tice in the flux flow region. With an external current, vortex 
in condensate with larger flux experiences stronger Lorentz 
force. On the other hand, the viscosity of vortex in differ- 
ent condensates is different due to the different size of normal 
core. As a result, vortices in some bands tend to move faster. 
For a small external current, the disparity of vortex motion 
can be balanced by the interband attraction between vortices 
in different condensates, and vortices in different condensates 
move with the same velocity, as shown in Fig. [T] However at 
a large current, the interband attraction may be insufficient to 
compensate the disparity of driving force and viscosity of vor- 
tex in different condensates and composite vortices are disso- 
ciated, i.e. vortices in different condensates move with differ- 
ent velocities. In the decoupled phase, the flux flow resistivity 
increases. The Shapiro steps are induced under an additional 
ac current when the oscillation of vortex lattice induced by 
the ac current is resonant with the vortex oscillations due to 
the periodic potential created by the relative motion of vortex 
lattice in different condensates. 

We use the London free energy functional for two-band su- 
perconductors with a Josephson-like interband coupling. The 
free energy density therefore can be written as 



+ (V x A) -y cos (0i -0 2 ) 

,— r- (1) 
where A M = ^(m^c 2 )/ (Am^e 2 ) is the London penetration 

depth for each condensate with superfluid density n M . A is vec- 
tor potential, is the electron mass in fi-th band and y is the 
interband Josephson coupling. The effective penetration depth 
for the two-band system is A~ 2 — 2^ =1 ^ n tne absence of 
external magnetic fields, the phase of different superconduct- 
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FIG. 1. (color online) (a) Schematic view of vortex lattice in two- 
band superconductors, (b): In the ground state, the normal core of 
vortex in one condensate with flux O2 (red circle) is locked with that 
in the other condensate with flux 0[ (green circle) to form a compos- 
ite vortex with one quantum flux ®o = ©2 + $>i- (c): With a small 
current, two vortex lattices are shifted with respect to each other but 
move with the same velocity and the composite vortex is deformed, 
(d) At a large current, the two vortex lattices are decoupled, move 
with different velocities, and the composite vortex is dissociated. 



ing condensates is locked, 8 X — 02 for 7 > and 6\ — 62 + n 
otherwise. The decoupling of phase due to the injection of 
current in superconducting wire both in equilibrium |[T2l and 
out of equilibrium [ 8 1 was discussed by Gurevich and Vinokur 
recently. 

Minimizing Tl with respect to A, we obtain the London 
equation 



(2) 



where O^, = A &0/A 2 is the fractional quantum flux and r^j = 
(x^j^n j) is the vortex coordinates for the vortex in the fi-th 
condensate. The vortex line is assumed to be straight. T L can 
be splitted into two contributions [20] Tl — T m + Tc with the 
magnetic coupling 



(3) 



T m is the same as that in single-band superconductors because 
there is only one gauge field A in superconductors. T m ac- 
counts for the magnetic coupling between vortices in different 
condensates. T c represents the coupling due to the phase dif- 
ference between condensates 

T c = [V (0! - 6 2 )] 2 - y cos (9, - 2 ) . (4) 

T c does not depend on A and accounts for the locking of 
two superconducting phases. For a fractional vortex where 



9\ changes by 27r around ro while 62 does not change, the 
self-energy per unit length is 

(5) 

where L is the linear size of the system and ^ is the coher- 
ence length. Ef v diverges at L — > 00 due to the neutral mode 
described by the term proportional to [V(0i - O2)] 2 in Eq. 
The Josephson contribution in Eq. <|3j renders the fractional 
vortex linearly divergent in L. Thus a fractional vortex is ther- 
modynamically unstable in bulk superconductors Q. For a 
composite vortex where 0„ changes by 27r around the same 
position, T c — and its self-energy is finite. 

To calculate intraband and interband interaction for vortex, 
one needs to know 9^. They can be obtained by minimizing 
Eq. Q with respect to 6\, 
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v z (e 1 -e 2 )-ysia(e 1 -e 2 ) = o, 



(6) 



subject to the boundary condition accounting for vortices 

Vx(V9 IJ ) = 2xY J 6(r-r lx>j ). (7) 



The interaction between vortices according to Eqs. |6]) and 
dVj is nonlinear, thus it is many-body interaction. The term 
y{Q\ - 62) is of the order 1/a 2 with a being the average dis- 
tance between vortices in the same condensate. For a strong 
field such that a <K Xj = -^/Oi<D 2 /(167r 3 ^i 2 |y|), the sine term 
becomes a 2 1 A 2 , times smaller than the gradient term thus can 
be neglected. In this case, T c reduces to the free energy for 
the XY model. For MgB 2 , y ~ 150 J/m 3 |[9]], it requires fields 
stronger than 4 T at temperature T = K. For V^Si ED and 
FeSei_ x [22|, the required field is smaller because the inter- 
band coupling y is much weaker. 

We then discuss the interband and intraband interaction be- 
tween vortices, neglecting the Josephson interband coupling 
term. Then both T m and T c are quadratic in B and 6^, therefore 
the interaction between vortices is pairwise. T m accounts for 
short-range interband and intraband repulsion between vor- 
tices with the same polarization. T c describes the plasma 
interaction in two dimensions and the interaction between 
vortices is long range. The term proportional to (V^) 2 ac- 
counts for intraband repulsion and the term proportional to 
-V9\V02 produces interband attraction between vortices in 
different condensates. This interband attraction outweighs the 
interband repulsion in T m - The intraband repulsion between 
two vortices in the same condensate separated by a distance 



— r u,i r u 



(8) 



8tt 2 A 2 "\ A I &7i 2 A 2 
and the interband attraction between two vortices in the dif- 
ferent condensates with a separation V[2jj = ry - r 2 j is 

o 1 o 2 



Vm er (r 12 jj) 



8n 2 A 2 



^) + ln(r 12 , y ) 



(9) 



3 



Equations ([8]l and |9]) are valid away from vortex cores. 

In the flux flow state, vortices in each condensate driven by 
the Lorentz force move and cause dissipation, resulting in vor- 
tex viscosity. The dissipation is due to the motion of normal 
core, thus the viscosity for vortices in each band depends on 
of the corresponding band. In the framework of the Bardeen- 
Stephen model, the viscosity is given by 77^ = <E> ( 2 /(27rc 2 ^ 2 ). 
We use the quasistatic approximation, i.e. the vortex structure 
in each condensate does not change in the dynamic region, 
and introduce overdamped dynamics for vortices 



<&1<$ 2 A 
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(10) 



where J is the external current. The effect of disorder be- 
comes less important in the flux flow region because the dis- 
order is averaged out by vortex motion, and lattice order- 
ing is improved 11231 l24l . In the lattice phase, the intraband 
vortex interaction vanishes due to symmetry. The interband 
vortex attraction can be written in the momentum space and 
we only take the contribution from the dominant wavevector 
G = (In/a, 0) for the vortex lattice moving along the x direc- 
tion. Here a is the lattice constant and we assume a square 
lattice. In the region 2nA/a » 1, the equation of motion for 
the center of mass of vortex lattice R^ in each band becomes 

i 2 d, (R 2 - Ri) = - (1 + rQ sin (R 2 - R{) + (G> 2 - rQ J, (11) 



d,Ri + i 2 d,R 2 = (1 + <DI,) J. 



(12) 



We have written Eq. ( 1 1 1 in term of the relative motion be- 
tween the vortex lattice in different bands, and the sine term 
accounts for the attraction between the two lattices. Here di- 
mensionless units are introduced: length is in unit of a/(2n), 
time is in unit of r\\a/(2nF d), current is in unit of cFd/$>\- 
Fd is the maximum attractive force between two lattices Fd = 
0!<I)2fl/(1287r 6 /l 4 ). <b' 2 = 2 /<I>i and = 772/771 . Similar 
equations for vortex motion in bilayer superconducting films 
was presented in Ref. ||25ll . 

At a small current, two vortex lattices in different 
band move as a whole with a velocity vi = v 2 = 

(l+ ?7 2 ) (l + <t> 2 ) J. The centers of mass of these two lat- 
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tices deviate with a separation sin 

At this stage, the composite vortex starts to deform. The max- 
imum attraction is reached at R 2 -Ri = n/2 or a/4 in real unit. 
At a threshold current 



h = \(J+n' 2 )(,%-ri 2 y 



(13) 



these two lattices decouple and move at different velocities, 
resulting in the dissociation of composite vortex. Their corre- 




FIG. 2. (color online) The dependence of velocity V\ , Vi and electric 
field E on the current J, obtained from Eqs. jTTJ ), ( | 1 2} and l |14[ >. 



sponding velocity is 
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(1 + V 2 )j _^J(o 2 -^) 2 /2_( 1+ ^) 2 



(14) 



The dependence of v M on J is displayed in Fig. [2] At a large 
current J » Jd, each lattice behaves independent with veloci- 
ties v\ - J and v 2 = J/t]' 2 - 

The decoupling of two lattices can be observed experimen- 
tally in transport measurements. The I-V characteristics can 
be derived from power balance condition 771V 2 + 772 v 2 = JEa 2 
with E the electric field. The I-V curve is shown Fig. |2j where 
the differential resistivity dE/dJ increases in the decoupled 
phase. Experimental observation of such increase may be 
challenging because the decoupling current usually is large, 
and the Larkin-Ovchinnikov instability of vortex lattice may 
be important|26]. The vortex core shrinks due to the electric 
field caused by vortex motion and the flux flow resistivity in- 
creases with current, which blurs the dissociation transition 
in the I-V curve. Nevertheless, the dissociation transition can 
be confirmed unambiguously by measurement of the Shapiro 
steps in the decoupled phase, as discussed below. 

In the decoupled phase, if one takes one lattice as reference, 
the other lattice experiences periodic potential induced by the 
reference lattice. When an ac current is added in superposition 
to the dc current, the oscillation of the moving lattice induced 
by the ac current may be in resonance with the oscillation due 
to the periodic potential of the reference lattice if the period of 
the ac current matches with the period of the potential. This 
results in Shapiro steps in I-V curves. The physics is the same 



as that in Josephson junctions and Eq. (Hi also describes the 
phase dynamics in overdamped Josephson junctions. With a 
current J — /dc + Re [/ ac exp[/(wf + <p)]\, the center of mass 
of each lattice is R^ = v^t + Re La^ exp[/(v2 - vi)f]J in the 



region |vi — v 2 \ = co » 1. From Eqs. (Hi and ( 12 1, we obtain 
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FIG. 3. (color online) The same as Fig. [2] but with an ac current in 
addition to a dc current. 
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vi) 771]. The dc current 



; exp(^)/Oi]/[(v 2 • 

is / dc = (<t>' 2 - ^ Re [rf 2 (v 2 - Vi) + (l + tf 2 ) (A 2 - AO/2]. 
When one changes /do <p adjusts correspondingly because vi - 
V2 is locked with the driving frequency to, and a Shapiro step 
is traced out. The height of the Shapiro step is 



(1+770(^-1) 
(v 2 -v 1 )(<5' 2 -77' 2 )' 



(15) 



The results for Shapiro step obtained by solving Eqs. ( 1 1 
and (12 1 with J — /d c + 1.2sin(o»f) numerically, are shown 
in Fig~[3] Shapiro steps appears when at = vi - v 2 . Here 
we only considered the dominant resonance. There are also 
Shapiro steps at na> = (v t - V2) with an integer n > 1 with 
a smaller height. In the presence of quenched disorder, the 
Shapiro steps occurs at nu> = v M (with the reduced units) due to 
the periodic passing of vortex lattice through the defects. J27] 
28 1 These steps can be distinguished from those induced by 
relative motion of two vortex lattices, because their resonance 
condition is different. 

The decoupling of composite vortex lattice depends on the 
two parameters < t>2/ < 3 , i an d 72/^1 • Generally the effect is 
present in all multiband superconductors. However for a small 
disparity between bands in ^ and A M , the current Jj is high. 
There are superconductors with large disparity among con- 
densates, such as the proposed liquid hydrogen superconduc- 
tors due to large mass difference between proton and electron. 
To estimate for MgB2, we take £1 = 13 nm, ^2 =51 
nm, A x = 47.8 nm and A 2 = 33.6 nm at T = K IfTTl . 
Pfx — 1CT 9 Q • m[29] and a — 40 nm corresponding to field 
at B « IT Then we obtain J4 — 5 x 10 9 A/m 2 , which is 
much smaller than the depairing current. The velocity of vor- 
tex lattice at dissociation transition is V\ — v 2 ~ 3 m/s, which 
is smaller than the typical Larkin-Ovchinnikov instability ve- 
locity for vortex lattice . QUI 



The dissociation of composite vortex lattice shares some 
similarity to that in multilayer superconductors PD and also 
in cuprate superconductors ll32l - l34l . In the latter case, vor- 
tices in different layers both carry <f>o flux and the dissoci- 
ation occurs in the real space, while for multiband super- 
conductors, vortices carry fractional flux and the dissocia- 
tion occurs in the band (momentum) space. The decoupling 
transition has been observed experimentally ||3T 1 43"4l and dis- 
cussed theoretically ||25l |35l |36l decades ago. The Shapiro 
steps are also observed in the decoupled phase in multilayer 
superconductors [37 1 . These observations in multilayer super- 
conductors corroborate the possible observation of the pre- 
dicted dissociation of composite vortex in multiband super- 
conductors. 

To summarize, we have predicted the dissociation of com- 
posite vortices in two-band superconductors in the flux flow 
region because of the disparity in vortex viscosity and flux of 
vortex in different condensates. At a small velocity, the two 
vortex lattices are shifted with respect to each other and the 
composite vortices are deformed. At a high velocity, the two 
vortex lattices move with different velocities resulting in the 
dissociation of composite vortices. In the decoupled phase, 
the flux flow resistivity increases. The Shapiro steps are in- 
duced when an ac current is applied in addition to the dc cur- 
rent. 
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